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Abstract 

Canonical forms for congruence and *congruence of square com- 
plex matrices were given by Horn and Sergeichuk in [Linear Algebra 
Appl. 389 (2004) 347-353], based on Sergeichuk's paper [Math. USSR, 
Izvestiya 31 (3) (1988) 481-501], which employed the theory of repre- 
sentations of quivers with involution. We use standard methods of 
matrix analysis to prove directly that these forms are canonical. Our 
proof provides explicit algorithms to compute all the blocks and pa- 
rameters in the canonical forms. We use these forms to derive canonical 
pairs for simultaneous congruence of pairs of complex symmetric and 
skew-symmetric matrices as well as canonical forms for simultaneous 
*congruence of pairs of complex Hermitian matrices. 

AMS classification: 15A21; 15A63 

Keywords: Canonical forms. Congruence, ^Congruence, Bilinear 
forms; Sesquilinear forms. Canonical pairs. 



1 Introduction 

Canonical matrices for congruence and *congruence over any field F of char- 
acteristic not 2 were established in |20l Theorem 3] up to classification of 
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Hermitian forms over finite extensions of F. Canonical forms for complex 
matrix congruence and *congruence are special cases of the canonical ma- 
trices in [20] since a classification of Hermitian forms over the complex field 
is known. Simpler versions of these canonical forms were given in 0, which 
relied on [2D], and hence on the theory of representations of quivers with 
involution on which the latter is based. 

In this paper, all matrices considered are complex. We use standard 
tools of matrix analysis to give a direct proof that the complex matrices 
given in [8] are canonical for congruence and *congruence. 

Let A and B be square complex matrices of the same size. We say that 
A and B are congruent if there is a nonsingular S such that S'^AS = B; 
they are ^congruent if there is a nonsingular S such that S*AS = B. We let 
S* := [sji] = 5^ denote the complex conjugate transpose of S = [sij] and 
write S-^ := (S"!)^ and S'* := {S''^)*. 

Define the n-by-n matrices 







i-ir 



-1 



1 1 

-1 -1 



1 1 



(ri = [1]), 



(1) 



1 

1 i 



(Ai = [1]), 



(2) 



and the n-by-n Jordan block with eigenvalue A 

"A 1 O" 



Jn(A) 



A 







(Ji(A) = [A]). 



The most important properties of these matrices for our purposes are that 
r„ is real and r~^r„ = r~*r„ is similar to Jri((— 1)"^^); A„ is symmetric 
and A~*A„ = A~^A„ is similar to Jn(l)- 
We also define the 2n-by-2n matrix 



In 



1 
fi 



(3) 
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the skew sum of Jn(/^) and /„. If /i 7^ 0, then ff2n(^) H2n{p) is similar to 
Jnilj) ® JnifJ-'^) and H2n{lj)^* H2n{lj) is similar to J„(/u) J„(/2~^). 

Sylvester's Inertia Theorem describes the *congruence canonical form of 
a complex Hermitian matrix. Our main goal is to give a direct proof of 
the following theorem, which generalizes Sylvester's theorem to all square 
complex matrices. 

Theorem 1 ([8, Section 2]). (a) Each square complex matrix is congruent 
to a direct sum, uniquely determined up to permutation of summands, of 
canonical matrices of the three types 



Type 


Jn(0) 


Type I 




Type II 


/i is determined up to replacement by 



(4) 



(h) Each square complex matrix is *congruent to a direct sum, uniquely 
determined up to permutation of summands, of canonical matrices of the 
three types 



Type 


Jn(0) 


Type I 


AA„, |A| = 1 


Type II 





Instead of An, one may use r„ or any other nonsingular n x n matrix Fn 
for which there exists a real On such that F~*Fn is similar to J„(e*^"). 

For *congruence canonical matrices of Type I, it is sometimes convenient 
to identify the unit-modulus canonical parameter A with the ray {tX : t > 0} 
or with the angle 9 such that A = e*^ and < ^ < 27r. If A occurs as a 
coefficient of exactly k blocks A„ in a *congruence canonical form, we say 
that it is a canonical angle (or ray) of order n with multiplicity k. 

Our proof of Theor em [1] provides explicit algorithms to compute the sizes 
and multiplicities of the canonical blocks Jn(0), r„, AA„, and H2n{lj) and 
their canonical parameters A and ^. 

It suffices to prove Theorem [T] only for nonsingular matrices because 
of the following lemma, which is a specialization to the complex field of a 
regularizing decomposition for square matrices over any field or skew field 
with an involution [9]. 
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Lemma 2. Each square complex matrix A is congruent (respectively, 
*congruent) to a direct sum of the form 



-B Jri (0) • • • © Jrp{0) with a nonsingular B. (6) 

This direct sum is uniquely determined up to permutation of its singular 
direct summands and replacement of B by any matrix that is congruent 
(respectively, *congruent) to it. 

The nonsingular direct summand in ([6]) is called the regular part of A; 
the singular blocks in ([6]) are its Type blocks. There is a simple algorithm 
to determine all of the direct summands in ([6]). If desired, this algorithm 
can be carried out using only unitary transformations [9]. 

In our development, it is convenient to use some basic properties of 
primary matrix functions. For a given square complex matrix A and a given 
complex valued function / that is analytic on a suitable open set containing 
the spectrum of A, the primary matrix function f{A) may be defined using 
a power series, an explicit formula involving the Jordan canonical form of 
A, or a contour integral. For our purposes, its most important property is 
that for each A, f{A) is a polynomial in A (the polynomial may depend on 
A, however), so f{A) commutes with any matrix that commutes with A. 
For a systematic exposition of the theory of primary matrix functions, see 
[6l Chapter 6]. 

Preceded by Weierstrass, Kronecker developed a comprehensive theory 
of equivalence of matrix pencils in the late nineteenth century, but a simi- 
larly complete theory of matrix congruence has been achieved only recently. 
Gabriel [3] reduced the problem of equivalence of bilinear forms to the prob- 
lem of equivalence of nonsingular bilinear forms. Riehm [16] reduced the 
problem of equivalence of nonsingular bilinear forms to the problem of equiv- 
alence of Hermitian forms. His reduction was improved and extended to 
sesquilinear forms in |17| . 

Using Riehm's reduction, Corbas and Williams [2] studied canonical 
forms for matrix congruence over an algebraically closed field with char- 
acteristic not 2. However, their proposed nonsingular canonical matrices 
are cumbersome and not canonical, e.g., their matrices 



A 



1 

1/2 



and B 



1 
2 



are actually congruent: BAB^ = B. For the singular case, they refer to 
the list of "singular blocks of known type" in [23^ p. 60]. These singular 
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blocks are canonical but cumbersome, and we are fortunate that they may 
be replaced by the set of singular Jordan blocks; see [20] or [9j . 

Any square complex matrix A can be represented uniquely as A = S + C, 
in which S is symmetric and C is skew-symmetric; it can also be represented 
uniquely as A = TC + ilC, in which both and /C are Hermitian. A si- 
multaneous congruence of S and C corresponds to a congruence of A and a 
simultaneous *congruence of Ti. and /C corresponds to a *congruence of A. 
Thus, if one has canonical forms for S and C under simultaneous congruence 
(often called canonical pairs), then one can obtain a canonical form for A 
under congruence as a consequence. Similarly, a canonical form for A under 
*congruence can be obtained if one has canonical pairs for two Hermitian 
matrices under simultaneous *congruence. Canonical pairs of both types 
may be found in Thompson's landmark paper [21] as well as in Lancaster 
and Rodman's recent reviews [2] and [13]. Thompson's canonical pairs 
were used to obtain canonical matrices for congruence over the real field by 
Lee and Weinberg [15j, who observed that "the complex case follows from 
Thompson's results just as easily." 

However, deriving canonical forms for complex congruence and * congru- 
ence from canonical pairs is like deriving the theory of conformal mappings in 
the real plane from properties of conjugate pairs of real harmonic functions. 
It can be done, but there are huge technical and conceptual advantages 
to working with complex analytic functions of a complex variable instead. 
Likewise, to derive congruence or *congruence canonical forms for a complex 
matrix A we advocate working directly with A rather than with its associ- 
ated pairs (5,C) or (7-^,/C). Our approach leads to three simple canonical 
block types for complex congruence rather than the six cumbersome block 
types found by Lee and Weinberg \15^, p. 208]. 

Of course, canonical pairs for (5,C) and (Ti,)C) follow from congruence 
and *congruence canonical forms for A. Define the following n-by-n matri- 
ces: 

1 

-1 ■•• 

■•• ■•• 1 

-10 



1 

1 





5 



(_l)n+r 




-1 



1 
-1 



Y 

1 ^ n 



10 
and a 2-parameter version of the matrix 

"0 

A„(a,6) : ^ 
a b 







1 
-1 
1 



-1)" 



a 

•■ b 





a,6 G C. 



Define the direct sum of two matrix pairs 

(^1, A2) e (^1,^2) := (^1 0^1,^2 B2) 
and the skew sum of two matrices 



[A\B] 



B 
A 



Matrix pairs (^1,^2) and {Bi,B2) are said to be simultaneously congru- 
ent (respectively, simultaneously *congruent) if there is a nonsingular matrix 
R such that Ai = R^BiR and A2 = R^ B2R (respectively, Ai = R*BiR 
and A2 = R*B2R). This transformation is a simultaneous congruence (re- 
spectively, a simultaneous *congruence) of the pair {Bi,B2) via R. 

The following theorem lists the canonical pairs that can occur and their 
associations with the congruence and *congruence canonical matrices listed 
in (jl]) and ([5]) of Theorem [H The parameters A and fi are as described in 
Theorem [H the parameters i/, c, a, and b in the canonical pairs are functions 
of A and fj,. 

Theorem 3. (a) Each pair {S,C) consisting of a symmetric complex 
matrix S and a skew- symmetric complex matrix C of the same size is 
simultaneously congruent to a direct sum of pairs, determined uniquely up 
to permutation of summands, of the following three types, each associated 
with the indicated congruence canonical matrix type for A = S -\- C: 
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Type 0: J„(0) 


(M„, Nn) 


Type I: r„ 


[Xn-,Yn) if n is odd, 
(Yn,Xn) if n is even 


Type II: i?2n(/") 
0/^/(-l)"+i 


([J„(/i + l)\J„(/i + lf], 

[J„(/i-l)\- J„(/i-in) 
and /i is determined up to replacement by ;U"^ 



(7) 



The Type II pair in ([7|) can be replaced by two alternative pairs 



Type II: if2n(/^) 

0//i/-l 
/i 7^ 1 if n is odd 


([Wn], [Jnil^)\-Ul^V]) 

u ^ if n is odd, v ^ ±1, 
u is determined up to replacement by —u 


Type II: H2n{-1) 
n is odd 


([J„(0)\J„(0)^], [In\-In]) 

n is odd 



(8) 



in which 

„='^. 

/i+ 1 

(b) Each pair (7Y,/C) of Hermitian matrices of the same size is simul- 
taneously *congruent to a direct sum of pairs, determined uniquely up to 
permutation of summands, of the following four types, each associated with 
the indicated congruence canonical matrix type for A = 7i + iJC: 



Type 0: J„(0) 


(M„, iNn) 


Type I: AA„, 
|A| = 1, AV -1 


±(A„(1,0), A„(c,l)) 
c G R 


Type I: AA„ 
A2 = -1 


±(A„(0,1), A„(1,0)) 


Type II: i^2n(/u) 
> 1 


([/„\/„], [Jn{a + ib)\Jn{a + ib)*]) 
a, 6 G M, a + bi^i, b> 



in which 

ImA 21m fj |/i|^ — 1 
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2 Congruence 

The cosquare of a nonsingular complex matrix A is A~'^A. If B is congruent 
to A, then A = S'^BS for some nonsingular S and hence 

^-^^ = {S^BS)-^{S^BS) = S-^B-^S~^S^BS = S-\B-^B)S, 

so congruent nonsingular matrices have similar cosquares. The following 
lemma establishes the converse assertion; an analogous statement for ar- 
bitrary systems of forms and linear mappings was given in [19j and [20\ 
Theorem 1 and § 2] . 

Lemma 4. Nonsingular complex matrices A and B are congruent if and 
only if their cosquares are similar. 

Proof. Let A'"^ A and B^^B be similar via S, so that 

A-^A = S-^B-^BS = {S-^B-^S-^) {S'^BS) = C-^C, (10) 

in which C := S^BS. It suffices to prove that C is congruent to A. Let 
M := CA-^ and deduce from (fTOl) that 



M = CA~^ = C^A-^ = {A-^C)^ and M'^ = A'^C. 

Thus, 

C = MA = AM^ 

and hence 

q{M)A = Aq{M^) = Aq{Mf 

for any polynomial q{t). The theory of primary matrix functions [H Section 
6.4] ensures that there is a polynomial p{t) such that p(M)^ = M, so p{M) 
is nonsingular and 

p{M)A = Ap{Mf. 

Thus, 

C = MA = p{MfA = p{M)Ap{M)^ 
so C is congruent to A via p{M). □ 

Proof of Theorem\^a). Let A be square and nonsingular. The Jordan 
Canonical Form of A^'^A has a very special structure: 

0(J^.(/x.)eJ^.(/z-i))e0j.,((-ir+^), o^^H^i-ir"-'; (11) 
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see \22\ Theorem 2.3.1] or [H Theorem 3.6]. Using (jlip . form the matrix 



i=l 



Since the cosquare 











J^(^)-^ J [JM 
is similar to Jmip-) © JmifJ-^^) and the cosquare 

T 



Jm(^) 

JM-^^ 



-L n -In 



(-1) 



n+1 



-1 



-1 



-1 



-1 



1 







(-1) 



n+1 



1 2 
1 





is similar to Jn((— 1)"^^), (jlip is the Jordan Canonical Form of B^^B. 
Hence A and B have similar cosquares. Lemma H] now ensures that A and 
B are congruent. Moreover, for any C that is congruent to A^ its cosquare 
C~'^C is similar to A~'^A, so it has the Jordan Canonical Form (jlip . which is 
uniquely determined up to permutation of summands. Hence B is uniquely 
determined up to permutation of its direct summands. □ 

Our proof of Theorem [T]^a) shows that the congruence canonical form of 
a square complex matrix A can be constructed as follows: 

1. Use the regularizing algorithm in [9l Section 2] to construct a regu- 
larizing decomposition ([6|) of A (if desired, one may use only unitary 
transformations in that algorithm). 

2. Let B be the regular part of A and determine the Jordan canonical 
form (jlip of its cosquare B~^B. 



3. Then 



H2m, (/"i) ® Tn, © Jri (0) © • • • © Jr^ (0) 



i=l 



is the congruence canonical form of A. 
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3 * Congruence 



The *cosquare of a nonsingular complex matrix A is A = A *A. If B is 
*congruent to A, then A = S*BS for some nonsingular S and hence 

A-*A = {S*BS)-*{S*BS) = S-^B-*S-*S*BS = S-^{B-*B)S, 

so *congruent nonsingular matrices have similar *cosquares. However, [—1] 
is the *cosquare of both [i] and [—i], which are not *congruent: there is 
no nonzero complex s such that —i = sis = Nevertheless, there is a 

useful analog of Lemma H] for *congruence. We denote the set of distinct 
eigenvalues of a square matrix X by dspecX. 

Lemma 5. Let A and B he nonsingular n-by-n complex matrices with sim- 
ilar *cosquares, that is, A~*A = S^^{B^*B)S for some nonsingular S. Let 
Bs ■= S*BS, let M := BsA~^ , and suppose M has k real negative eigenval- 
ues, counted according to their algebraic multiplicities (0 < k <n). Then: 

(a) M is similar to a real matrix. 

(b) There are square complex matrices D_ and of size k and n — k, 
respectively, such that A is *congruent to (—D^) ©D+ and B is *congruent 
to D^®D+. 

Proof. We have 

A-*A = S-\B-*B)S = {S-^B-*S-*) {S*BS) = Bs*Bs, 
from which it follows that 

M = BsA-^ = B*sA-* = {A-'^BsY and M* = A-'^Bs, 
and hence 

Bs = MA = AM*. (12) 

Thus, M* = A~^MA, so M is similar to a real matrix [5l Theorem 4.1.7]. Its 
Jordan blocks with nonreal eigenvalues occur in conjugate pairs, so there is 
a nonsingular T such that TMT~^ = M_ ©M_|_, in which the k-hy-k matrix 
M_ is either absent or has only real negative eigenvalues; has no negative 
eigenvalues and is similar to a real matrix if it is present. Moreover, if we 
partition TAT* = [Aijljj^-^ conformally to M_ © M+, Sylvester's Theorem 
on linear matrix equations [6, Theorem 4.4.6] ensures that 

TAT* = Au(BA22 (13) 
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since the equalities dspec M_ n dspec M4 
dspec M+ = dspec Ml imply that 



0, dspec M_ = dspec iVP, and 



TBsT* = {TMT~^) (TAT*) = (TAT*) {T-*M*T*) 
= (M_ e M+) (TAT*) = (TAT*) (M* Ml 



M^Aii M^Ai2 
M+A21 M+A22 



AiiMl A12M* 
A21MI A22M* 



Thus, 



so 



and 



= M_Aii e M+A22 = AuMl A22MI. 

{-M^)]Au = Aui-Ml) and M+A22 = A22M; 
qi{-M_)An = Anqi{-Ml) = Auqi{-M^y 



q2{M+)A22 = A22q2{MX) = A22q2{M+y 

for any polynomials qi{t) and q2{t) with real coefficients. Neither — M_ nor 
-M+ has any negative eigenvalues and each is similar to a real matrix, so O 
Theorem 2(c)] ensures that there are polynomials g{t) and h[t) with real 
coefficients such that g{-M_f = -M^ and h{M+f = Af+. It follows that 
B is * congruent to 



TBsT* 



-M^Au) e M+A22 = -giM^fAii h{M+fA22 



-g{-M+)Ang{-M+)* h{M+)A22h{M+r 



in which D_ = -g{-M+)Aiig{-M+)* and D+ = h{M+)A22h{M+)*; A is 
*congruent to -D+ by ([13]) and i? is *congruent to □ 

The result cited from [71 Theorem 2(c)] gives sufficient conditions for 
f{X) = y to have a real solution X for a given real y. The key conditions 
are that / is analytic and one-to-one on a domain that is symmetric with 
respect to the real axis and is typically real, that is, f{z) = f{z) on the 
range of /. Under these conditions there is a solution X that is a polynomial 
in Y with real coefficients. In the case at hand, f{z) = on the open right 
half-plane; this special case appears in [10^ p. 545] and was employed in [211 
p. 356] and [14', Lemma 7.2] to study canonical pairs of Hermitian matrices. 

Two noteworthy special cases of Lemma [5] occur when M either has only 
positive eigenvalues (k = n) or only negative eigenvalues {k = 0). In the 
former case, A is *congruent to B; in the latter case, A is *congruent to 
-B. 
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Proof of Theorem\^h): Existence. Let A be nonsingular and let A = A^*A 
denote its *cosquare. Because A~* = {A~*A)~* = AA"* = AAA~^, for 
each eigenvalue A of ^ and each /c = l,2,..., Jfc(A) and Jfc(A^^) have equal 
multiplicities in the Jordan Canonical Form of A. Since A = A~^ whenever 
|A| = 1, this pairing is trivial for any eigenvalue of A that has modulus one; 
it is nontrivial for eigenvalues A whose modulus is greater than one. 
Let 



> 1, < (t>j < 2tt (14) 



i=l 



be the Jordan Canonical Form of A and use it to construct the matrix 

B := 0F2„,(m.) e 0e'<^^/2A„^,. (15) 



i=l 



Since the *cosquare 









Jm{^J) 







JM-* 



is similar to JmifJ-) © ^m(/U ) and the *cosquare 



1 2i 
1 





2i 
1 



(16) 



is similar to Jn(l), (|14p is also the Jordan Canonical Form of B~*B. Hence 
A and B have similar *cosquares, whose common Jordan Canonical Form 
contains 2p + q Jordan blocks {p block pairs of the form JmifJ-) © Jmifi'^) 
with |//| > 1 and q blocks of the form AA„ with |A| = 1). 

Lemma [5] ensures that there are matrices and such that A is 
*congruent to (—!?_) © and B is *congruent to D_ © D^. If D- is 
absent, then A is *congruent to B; if is absent, then A is *congruent to 
—B. In both of these cases A is *congruent to a direct sum of the form 



0ei^2™,(/i^)©0<Jie*'^^/2^„^ 



ei,6j e {-1,+1}. 



(17) 



i=l 



If both direct summands and are present, then their sizes are less 
than the size of A. Reasoning by induction, we may assume that each of 
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Z)_ and is *congruent to a direct sum of the form (|17p . Then A is 
*congruent to a direct sum of the form (jl7p as well. We may take all £i = 1 
in p!7|) since each H2m{P') is *congruent to —H2mifJ')- 



' Im ■ 




"0 ■ 








" 


_ -/^ 




Jm(^) 




. -Im _ 




Jm{^J) 



□ 

We have demonstrated that a nonsingular A is *congruent to a direct 
sum of Type I and Type II blocks 

p <? 

0^2m,(/^i) e 0(5,e^<^^/2A,^, 5,- G {-1, +1}, l^il > 1, < (l>j < 27r, 
1=1 j=i 

(18) 

in which the sizes 2mj and parameters /Xj of the Type II blocks as well as 
the sizes rij and squared parameters (Jje*'^^/^)^ = e**^^ of the Type I blocks 
are uniquely determined by A. Our reduction algorithm using Lemma [5] 
determines a set of signs {6j} that gives the desired *congruence of A to 
(jlSp . but we must show that no other choice of signs is possible: in the set 
of Type I blocks in (jlSp with equal sizes nj and equal coefficients e**^^/^ the 
number of blocks with sign equal to +1 (and hence also the number of blocks 
with signs equal to —1) is uniquely determined by A. 

Proof of Theorem\^h): Uniqueness. Let each of A and be a direct sum 
of Type I and Type II blocks and suppose that A and B are *congruent. 
We have shown that A and B have the form 

p q p q 

^ = ^2m. (/Xi) e Xj^n, , ^ = H2m, (fn) A„ . , 

i=l j=l 1=1 j=l 

in which all Hj € {— 1, +1}, all > 1, and all |Aj| = 1. Our goal is to 
prove that each of these direct sums may be obtained from the other by 
a permutation of summands. We may rearrange the summands to present 
A = Ai®A2 and B = Bi® B2, in which 

k k 
Ai = ^^Ar.A„j,, Bi = ^^KrXr^rirJ 
r=l r=l 
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and 



Xl / Xj for all . 



k + 1, . . . ,q. Let S*AS = B and partition 



S = [Sij\^ conformally with Ai A2. Since the *cosquares of A and B 
are similar via S, we have S {B~*B) = {A~*A) S and hence 



Su {B^ *Bi) S12 {B2 *B2) 
S21 {Bi *Bi) S22 {B2 *B2) 



{A^*Ai)Su {A^*Ai)Si2 
(^2 *^2) ^21 {A2*A2) S22 



But 



and 



dspec *B2) n dspec {A^ *Ai 







dspec {B^ *Bi) n dspec {A^ M2) = 0, 

so Sylvester's Theorem on linear matrix equations ensures that S = Su © 
522- Thus, Ai is *congruent to Bi via and hence it suffices to consider 
the case A = Ai and B = Bi. Moreover, dividing both A and B by Ai it 
suffices to consider a pair of *congruent matrices of the form 

k k 

A = 0e,A„,„5 = 05,A„^, € {-l,+l}. (19) 



r=l 



r=l 



We may assume that the summands in (I19p are arranged so that 1 < 
ni < • ■ ■ < rik. Define 



Let 



denote 



and let 



in which 



N := {ni, . . . ,nk) and 

Jn := Jni(O)©---© Jnfe(O) 

Aat := A„i ©••■ ©A„j^ and Pat : = 



ni H h nfc. 



(20) 
(21) 



N 



Pr, 



P 



ni 



P 



is the n-by-n reversal matrix. 

Since Pn^n = I + iJn and Jn is nilpotent. 



V 



N 



A-*An 



{PnA*^)-^ PnAn = (/ - iJN)-\l + iJTv) 



(/ + iJN + i'^JN + i^JN + ■■■){! + Un) 



I + 2iJN + 2i'^jl + 2i^Jif + 



,■3 t3 
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is a polynomial in Jtv- Moreover, J^r = i{I — 'Dn){I +'Dn)~'^ and (/+2?Ar)~^ 
is a polynomial in Pat, so Jat is a polynomial in Pat. 

Let C be |A^|-by-|A^| and partition C 
each block Cij is ni-hy-rij. Then C commutes with J^r if and only if each 



[^ij]i,j=i conformally to Jj\r, so 



block Cij has the form 



c 



.(2) 



C 



if i < j, 



.(2) 



Cii 



.(2) 



(22) 



if i > j; 







see [U Section VIII, §2] or [6l Lemma 4.4.11]. Each diagonal block Ca is 
upper Toeplitz; each block Cij with i < j has an upper Toeplitz submatrix 
that is preceded by a zero block; each block Cij with i > j has an upper 
Toeplitz submatrix with a zero block below it. An |A^|-by-|A^| matrix whose 
blocks have the form (j22p is said to be N -upper Toeplitz. 

If C is A^-upper Toeplitz, then so is PnC*Pn'-, the upper Toeplitz subma- 
trix in its i, J block is the complex conjugate of the upper Toeplitz submatrix 
in the j, i block of C. The matrices Pat, Pn^, and PnB (see (fT9]l ) are all 
A^-upper Toeplitz. 

Let C = [Cjj]^j^]^6e a given A-upper Toeplitz matrix. Consider the 
mapping C ^ C_ that takes C into the k x k matrix whose i,j entry is Cij 
if rii = rij and is otherwise; see ()22p . Partition the entries of N (already 
nondecreasingly ordered) into groups of equal entries 



1 < ni 



ni < ni+i 



nt < 



and observe that C_ is structurally block diagonal: 

c = Ci e C2 e Cs e • • • 



(23) 
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The sizes of the direct summands of C_ are the multiphcities of the entries 
of N, that is, 



Ci :-- 



Cll 



Crl 



Clr 



Co 



Cr+l,r+l 



Ce,r+1 



Cr+ll 



(24) 



In addition, for any A^-upper ToepUtz matrix D we have CD = C_ ■ D.- If 
C is nonsingular, we have / = CC~^ = C_ ■ C~^ , which imphes that C_ is 
nonsingular. A computation reveals that PnC*Pn = iC^* ■ 



Finally, consider the *congruent matrices A and B in (jl9p . which satisfy 
A~*A = B-*B = Vn. We have PnA = diag(ei, . . . , e^) and PnB = 
diag((5i, ...,6k). Let S*AS = B, so B-*B = S'^ {A-*A) S and hence 

S{B-*B) = SVn = VnS = {A~*A) S. 

Since S commutes with TDn and Jn is a polynomial in T>m, S commutes 
with Jat and hence is A^-upper Toeplitz. Moreover, 

PnB = PnS*AS = {PnS*Pn) {PnA) S 

and each of the matrices PnB, PnA, S, and PnS*Pn is A^-upper Toeplitz. 
Therefore, 

PnB = ( PnS*Pn ) (PnA ) S, 

that is, 

diag(5i, ...,6k) = S_* diag(ei, . . . , ek)S. 
But 5 is nonsingular and block diagonal: 

S = Si® 82^83® ■■■ , 

in which the respective nonsingular blocks 81,82,.-. are the same size 
as the respective blocks Ci,C2,... in (fM|) . Thus, diag{6i, . . . ,Sr) = 
S'J'diag(ei,...,e^)S'i, diag{5r+i, . . . ,Se) = 8^diag{£r+i, ... ,81)82, etc. 
Sylvester's Inertia Theorem ensures that diag((5i, . . . , 5^) can be ob- 
tained from diag(ei, . . . , by a permutation of its diagonal entries, 
diag((5r+i, . . . ,61) can be obtained from diag(er+i5 • • • , e^) by a permutation 
of its diagonal entries, etc. Therefore, each of the direct sums ([T9]) can be 
obtained from the other by a permutation of summands. □ 

The argument that we have just made also clarifies the final assertion in 
Theorem[T] A *cosquare F~*Fn is similar to Jn(A) with |A| = 1 if and only 
if it is not decomposable into a nontrivial direct sum under similarity. 
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4 An alternative algorithm for *congruence 



Although we can now determine the *congruence canonical form of a non- 
singular complex matrix A, in practice it is useful to have an alternative 
algorithm. 

Let A = A~*A, let /Ui,...,/!^ be the distinct eigenvalues of A with 
modulus greater than one, and let Ai, . . . , be the distinct eigenvalues of 
A with modulus one. Let S be any nonsingular matrix such that 



in which dspecCj = {^i,fL~ } for i = 1, . . . , r and dspecCr-+i = {Aj} for 
i = 1, . . . ,s. One way to achieve this decomposition is to group together 
blocks from the Jordan Canonical Form of A, but other strategies may be 
employed. Partition S*AS = [Aij]l~j^-^ conformally to the direct sum in ([25]) . 
The argument in the proof of uniqueness in Section [3] shows that S*AS is 
block diagonal: 



in which each An is the same size as Cj, dspecA~-*Aii = {fii,fL^ } for 
i = 1, . . . , r, and dspec A'^* An = {Xi^r} for i = r + 1, . . . , r + s. 

The Type II *congruence blocks are now easy to determine: to each 
pair of Jordan blocks Jm(^j) © JmiP-^^) of A~* An corresponds one Type II 
*congruence block -ff2m(/^i) of A. 

Now consider each diagonal block Ay+j^r+j in turn. Its *cosquare has a 
single eigenvalue \j = e'^'^^ , < 4>j < 2tt. Let the Jordan Canonical Form 
of the *cosquare of e~^'^i^'^Ar-^-j^r+j be / + Jn, in which := (ni, . . . ,nk) 
and 1 < ni < ••• < Uk- That *cosquare is similar to Pat := ^Ji*^N] 
let S be nonsingular and such that SV]\fS~^ = e~^'^^ A'J^^ -j^^Ajj^r-j+r- For 
notational convenience, normalize and set A := e~^^i ^'^ S* Aj^r,j+rS ■ Then 



— 1 



(25) 



S AS — All ■ ■ ■ ® Aj-r ^r+l,r+l © ' ' ' ^s+l,s+l 



(26) 



A-*A = e-'^^S-^A^ 



S = Vn = Aj^*A]sf, 



which implies that A = A*A 



■N 



*An and hence 



B := A],' A = A],\A*A],*An) = AjKAJ^'AYAn 
= A]^^B*An = (A^i3A^*) An = V'^^EVn. 



(27) 



Thus, Vn commutes with B, so B is A- upper Toeplitz. 
Invoking the identity 

B* = A^BA^ , 
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already employed in ([27|) . compute 

PnB*Pn = Pn (AatBA^I) Pn = (PnAn) BA-^P^. 

Since P/v^Af = / + iJn commutes with any A^-upper Toeplitz matrix, it 
follows that 

PnB*Pn = B (PnAn) A^iPjv = BP^ = B. (28) 

Let X = [Xij] be a given A^-upper Toeplitz matrix, partitioned as in 
(j22p . The N -block star of X is the A^-upper Toeplitz matrix 

:= PnX*Pn. 

The upper Toeplitz submatrix of each i,j block of X^ is the complex conju- 
gate of the upper Toeplitz submatrix of Xji. We say that X is N-Hermitian 
if X^ = X, in which case the upper Toeplitz submatrices of each pair 
of blocks Xij and Xji are complex conjugates. If = (1, 1, . . . , 1) then 
X^ — X* and X is A^-Hermitian if and only if it is Hermitian. 

The identity (j28p asserts that B is A^-Hermitian. 

If X and Y are A^-upper Toeplitz, one checks that 

(XY)^ = Y^X^. 

We say that A^-upper Toeplitz matrices X and Y are ^congruent {N -block 
star congruent) if there exists a nonsingular A^-upper Toeplitz matrix S 
such that S^BS = C; ^^ongruence is an equivalence relation on the set of 
A^-upper Toeplitz matrices. 

Since B is A^-upper Toeplitz and A^-Hermitian, for any A^-upper Toeplitz 
matrix S we have 

S*AS = S* (AnB) S = {PnS^Pn) AnBS = PnS^{PnAn)BS 

= Pn{PnAn)S^BS = An (S^BS) . (29) 

If we can find a nonsingular A^-upper Toeplitz S such that 

S^BS = ei/„, e • • • e Eklrik with Si = ±1, 

it follows from (I29p that A is *congruent to 

eiAn, ® ■ ■ ■ e EkAn^. (30) 

Theorem [Dj^b) ensures that for each n = 1,2,... there is a unique set of signs 
associated with the blocks An of size n in (j30p . The following generalization 
of Sylvester's Inertia Theorem provides a way to construct these signs. 
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Lemma 6. Let C be nonsingular, N-upper Toeplitz, and N -Hermitian. 
Then there is a nonsingular N-upper Toeplitz matrix S such that 

S^S 



£klnk 



each Ei £ {—1, !}• 



(31) 



Proof. Since C is nonsingular, C_ and hence all of the direct summands in 
(|23p are nonsingular as well. 

Step 1. If cii 7^ 0, proceed to Step 2. If cn = 0, then cij ^ for some 
j < r since Ci is nonsingular. Let Sg = [Sij]jlj^i be the A^-upper Toeplitz 
matrix in which the diagonal blocks are identity matrices and all the other 
blocks are zero except for Sji := e^^In^ 
1, 1 entry of ^CSg is 



e*^/„, for a real 6 to be determined. The 



e'^cij + e' 



{cji = cij and cjj 
nonzero. 



'^cij + e'^e' 
since ■ 



■ie 



-2Re{e'''cij) + Cjj (32) 
C). Choose any 9 for which ([32]) is 



Step 2. We may now assume that cn is a nonzero real number. Let 
a := [ciil^^/^, so a^cii = ±1. Define the real A^-upper Toeplitz matrix 
S = al and form S^CS 



SCS, whose 1, 1 entry is ±1. 
Step 3. We may now assume that dspecCn = {cn} 



{±1}. Then 



dspec(ciiCj^j^^) = {1}, so there is a polynomial p{t) with real coefficients 

is upper Toeplitz 
Uni- Define the 



such that p{Cr 



"2 



nir = ciiCfi^ B Theorem 6.4.14], p{Cuj 
and commutes with Cn, and p{Cn)Cup{Cn) = cnl = 
real A^-upper Toeplitz matrix S = p(Cii) /, 
whose 1, 1 block is ±/„^. 

Step 4. We may now assume that Cn 
Toeplitz matrix 



ni ■ 



. and form S^CS, 
Define the A'^-upper 



S 
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and form S^CS; its 1, j and j, 1 blocks are zero for all j = 2, . . . , n^. 

The preceding four steps reduce C by ^congruence to the form ±/„^ ©C". 
Now reduce C in the same way and continue. After k iterations of this 
process we obtain a real diagonal matrix 



' ^klrik, G {-1, 1} 



that is t^congruent to the original matrix C. 



(33) 
□ 
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Thus, to determine the *congruence canonical form of a given square 
complex matrix A, one may proceed as follows: 

1. Apply the regularization algorithm [U] to determine the singular *con- 
gruence blocks and the regular part. This reduces the problem to 
consideration of a nonsingular A. 

2. Let S be any nonsingular matrix that gives a similarity between the 
*cosquare of A and the direct sum in (|25p and calculate S*AS, which 
has the block diagonal form (j26p . Consider each of these diagonal 
blocks in turn. 

3. Determine the Type II *congruence blocks by examining the Jordan 
Canonical Forms of the diagonal blocks whose *cosquare has a two- 
point spectrum: one Type II block H2mifJ'i) corresponds to each pair 
JmiP'i) © Jm{U^^) in the Jordan Canonical Form of A~* An. 

4. For each diagonal block Ajj whose *cosquare has a one-point spectrum 
(suppose it is e^'^), consider Ajj = e~'^'^^'^ Ajj. Find an S such that 
A~*Ajj = 5~^(A^^Aiv)5 and consider B = ll^S*AjjS, which is 
A'^-upper Toeplitz and A-Hermitian. 

5. Use Lemma [6] (or some other means) to reduce B by ^^ongruence 
to a diagonal form ([33]) . Then Ajj is *congruent to a direct sum of 
the form (I30p and the diagonal block Ajj corresponds to a direct sum 
eie*'^/2^„i © • • • © efce''^/2A„^ of Type I blocks. 

5 Some special *congruences 

A square complex matrix A is diagonalizable by *congruence if and only 
if the *congruence canonical form of A contains only 1-by-l blocks (which 
can only be Type blocks Ji(0) = [0] and Type I blocks AAi = [A] with 
|A| = 1). The Type *congruence canonical blocks for A are all 1-by-l if 
and only if A and A* have the same null space; an equivalent condition is 
that there is a unitary U such that 

A = U*{B (B Ok)U and B is nonsingular. (34) 

The *congruence class of the regular part B is uniquely determined, so the 
similarity class of its *cosquare B~*B is also uniquely determined. There 
are no Type II blocks for A and its Type I blocks are all 1-by-l if and only if 
the *cosquare of its regular part is diagonalizable and has only eigenvalues 
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with unit modulus. Thus, A is diagonalizable by *congruence if and only if 
both of the following conditions are satisfied: (a) A and A* have the same 
null space, and (b ) the *cosquare of the regular part of A is diagonalizable 
and all its eigenvalues have unit modulus. 

Let A be nonsingular and suppose that its *cosquare A is diagonalizable 
and all its eigenvalues have unit modulus. Let S be any nonsingular matrix 
that diagonalizes A, and consider the forms that ()25p and ()26p take in this 
case: 

A-*A = S{Ci®---®Cq)S-^ 

and 

S*AS = Ei®---(BEq, 

in which Cj = e^*^^/ for each j = 1, ...,q, < 9i < ■ ■ ■ < 6q < n, each 
Ej is the same size as Cj, and each EJ*Ej = e^*^^/. We have e~^^^Ej = 
giOj ^* — (e~*^j£'j)*, so each matrix e~^^^ Ej is Hermitian. Sylvester's Inertia 

Theorem ensures that e~^^^ Ej is *congruent to -^„+ffi(~-^„-) for nonnegative 

j j 

integers n'j and nJ that are determined uniquely by (the *congruence class 
of) e~^'^^Ej. It follows that each Ej is *congruent to 

e''^I^+ e {-e^'^I^-) = e*'^7„+ (e'^'^+^h^-), 

j j j 3 

SO A is *congruent to the uniquely determined canonical form 




< ^1 < • • • < 0g < vr. The angles 9j for which > 1 together with the 
angles 9j + TT for which nJ > 1 (all 0j € [0, vr), j = 1, . . . , g) are the canonical 
angles of order one of A; the corresponding integers and nJ are their 
respective multiplicities. 

We have just described how to determine the signs ej that occur in (j30p 
when = (1, . . . , 1); in this special case, every |A^|-by-|A^| matrix is A^-upper 
Toeplitz, and A^-Hermitian matrices are just ordinary Hermitian matrices. 
Two square complex matrices of the same size that are diagonalizable by 
*congruence are *congruent if and only if they have the same canonical 
angles of order one with the same multiplicities. 

This observation is a special case of a more general fact: over the re- 
als or complexes, each system of forms and linear mappings decomposes 
uniquely into a direct sum of indecomposables, up to isomorphism of sum- 
mands [20^ Theorem 2]. This special case was rediscovered in [12], which 
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established uniqueness of the canonical angles and their multiplicities but 
did not determine them (the signs £j remained ambiguous) except in special 
circumstances, e.g., if the field of values of e^'^A lies in the open right half 
plane for some (p G [0, 27r). Using (j34p to introduce generalized inverses and 
a natural generalized *cosquare, [T5] later gave an alternative approach that 
fully determined the canonical angles and their multiplicities for a complex 
matrix that is diagonalizable by *congruence. 

If A is normal, it is unitarily *congruent to A©Ofc, in which A is a diagonal 
matrix whose diagonal entries are the nonzero eigenvalues of A (including 
multiplicities). The canonical angles (of order one; no higher orders occur) 
of a normal matrix are just the principal values of the arguments of its 
nonzero eigenvalues; the multiplicity of a canonical angle is the number of 
eigenvalues on the open ray that it determines. Thus, two normal complex 
matrices of the same size are *congruent if and only if they have the same 
number of eigenvalues on each open ray from the origin. This special case 
was completely analyzed in [llj . If A is Hermitian, of course, its nonzero 
eigenvalues are on only two open rays from the origin: the positive half- 
line and the negative half-line. For Hermitian matrices, the criterion for 
*congruence of normal matrices is just Sylvester's Inertia Theorem. 

Finally, suppose that the *congruence canonical form of a nonsingular 
n-by-n complex matrix A has only one block; the singular case is analyzed 
in [9]. If n is odd, that block must be a Type I block AA„ (with |A| = 1). 
If n = 2m, it can be either a Type I block or a Type II block H2m{l^) (with 

> 1). Of course, the latter case occurs if and only if the *cosquare of A 
is similar to JmitJ-) © Jm{ji~^) and > 1. 

Suppose the *cosquare of A is similar to Jn{^) with [A[ = 1. Then the 
*congruence canonical form of A is eAA„ and e G {1,-1} can be determined 
as follows. Let A = X'^A, let S be such that A^^A^ = S-^{A-*A)S, 
and let M := S*ASA^^. Lemma [5]^a) tells us that M is similar to a real 
matrix. Because A is indecomposable under *congruence, in Lemma ^Mjo) 
either k = n or k = 0: 

{ — 1 if all the eigenvalues of M are negative 
(35) 
1 if no eigenvalue of M is negative. 

Alternatively, we can employ the algorithm described in Section HI Ex- 
amine B := A~^S*AS, which must be nonsingular, upper Toeplitz, and 



22 



real, so 



B 



6n bf^ 







hi 



^(2) 
"11 

foil 



The reduction described in Lemma [6] is trivial in this case, and it tells us 
that e is the sign of bu. 



Example 7. Consider 
whose *cosquare 



A 



A-*A 



1 2 
1 



1 2 
-2 -3 



is not diagonal, has —1 as a double eigenvalue, and hence is similar to 
J2(— 1). Thus, A is *congruent to eiA.2 with e = ±1. Let A = —iA and 
verify that 



1 2i 
1 



A2 *A2 = S-^{A-*A)S for 5 



-1 
1 i 



Then both eigenvalues of M = S*AS/S.2^ = —I2 are negative, so (j35]) ensures 
that e = — 1 and A is *congruent to —iA2. Alternatively, B = A~^S*AS = 
—I2, so the sign of 611 is negative and e = —1. 



Example 8. Consider 
whose *cosquare 



A 



A-*A 



1 
-1 1 



-1 2 
-1 



(36) 



is similar to J2(— 1). Thus, A is *congruent to £iA2. Let A = —iA and 
verify that 



1 2i 
1 



A^*A2 = S-^{A-*A)S for 5 



1 

-i 



Then both eigenvalues of M = S*ASA2 ^ = —I2 are negative, so (|35]) ensures 
that e = — 1 and j4 is *congruent to — 'iA2. 
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Example 9. Suppose |A| = 1 but / —1. Let a denote the real part of A 
and consider 



A 






A/a ■ 


_ A 


" 1 " 


A/a 


i 


a 


1 aAi 



whose *cosquare 



aAz 1 




" 


1 


= A2 


" 1 2o?i ' 


1 




1 


aAi 


1 



A-*A = y 



is similar to J2(A^). Thus, A is *congruent to eX/S.2. Let A = \ ^A and 
verify that 



1 2i 
1 



■ I\-^*/\2 = S-'^{A-*A)S for 5:= 
Then both eigenvalues of 



l/a^ 
l/a3 \i/a^ 



-i 1 
1 



1 

l/a2 



l/a^ 



have the same sign as a. Thus, ([35]) ensures that A is *congruent to AA2 if 
Re A > and to -AA2 if Re A < 0. 



6 Canonical pairs 

We now explain how to use the canonical matrices in Theorem [1] to obtain 
the canonical pairs described in Theorem [3l 

Proof of Theorem\3(a) . Each square matrix A can be expressed uniquely as 
the sum of a symmetric and a skew-symmetric matrix: 

A = S{A)+C{A), S{A):=\{A + A^) , C{A) := \ {A - A^) . 

Since S{R^ AR) = R^S{A)R and C{R^AR) = R^C{A)R, any congruence 
that reduces A to a direct sum 

R^AR = Bi®---®Bk 

gives a simultaneous congruence of {S{A),C{A)) that reduces it to a direct 
sum of pairs 

{S{Bi),C{Bi)) e • • • e {S{Bk),C{Bk)) . 
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Theorem HJa) ensures that A is congruent to a direct sum of blocks of the 
three types 

2J„(0), r„, 2HM, (38) 

in which 7^ /u / (—1)"'+^ and ^ is determined up to replacement by fi~^, 
and that such a decomposition is unique up to permutation of the direct 
summands. 

Computing the symmetric and skew-symmetric parts of the blocks (|38p 
produces the indicated Type 0, Type I, and Type II canonical pairs in ([7]). 

It remains to prove that the two alternative pairs in ([8]) may be used 
instead of the Type II pair 

([J„(/i + l)\J„(/i + 1)^], [J„(^ - 1)\ - Jn(^ - 1)^]) . (39) 

First suppose that = —1, so n is odd (since ^ / (—1)"+-'^) and we have 
the Type II pair 

{[JniO)\Jn{Of], [Jn(-2)\ - J„(-2)^]) . 

A simultaneous congruence of this pair via 

" In 

_ J„(-2)-^ 

transforms it to the pair 

(\ Jn{0fJni-2y 

\[ J„(-2)-ij„(0) 

Since Jn(— 2)~^ J„(0) is similar to Jn(0), there is a nonsingular matrix S 
such that S~^Jn{0)S = J„(— 2)~^ J„(0). Then a simultaneous congruence of 
(|40]1 via 

\ s-^ " 

5^ 

transforms it to the second of the two alternative pairs in ([8|). 

Now suppose that /U 7^ — 1. Let S he a nonsingular matrix such that 

S-'jn{fl-l)Jn{fl+l)-'S=Jn{u), ■= ■ (41) 

/U + 1 

A simultaneous congruence of the pair (j39p via 

"J„(/i + l)-i5 " 
5-^ 





-/„ 


1 


In 



(40) 
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transforms it to the first alternative pair in The definition ()4ip ensures 
that u ^ 1; u ^ —1 since fi ^ 0; and 7^ if n is odd since fi 7^ (— 1)"^'^. 
Because /i is determined up to replacement hy fi~^, v is determined up to 
replacement by 

//-i - 1 _ 1 - /i _ 

JF^TT ~ i + fi ~ 

□ 

Proof of Theorem\3(h) . Each square complex matrix has a Cartesian de- 
composition 

A = n{A)+ilC{A), n{A):=\{A + A*), 1C{A) := '^{-A + A*) {A2) 

in which both ^{{A) and fC{A) are Hermitian. Moreover, if Ti' and /C' are 
Hermitian matrices such that A = fi'+iJC' , then Tl' = ^{{A) and /C' = IC{A). 
Since H{R*AR) = R*n{A)R and IC{R*AR) = R*IC{A)R, any *congruence 
that reduces A to a direct sum 

R*AR = Bi®---®Bk 

gives a simultaneous *congruence of (7^(j4), 1C{A)) that reduces it to a direct 
sum of pairs 

{n{Bi)MBi)) e • • • e {n{Bk)MBk)) ■ 

Theorem WL^) ensures that A is congruent to a direct sum of blocks of the 
three types 

2J„(0), AA„, and H2n{fJ'), in which |A| = 1 and \fi\ > 1 (43) 

and that such a decomposition is unique up to permutation of the direct 
summands. The Cartesian decomposition of 2J„(0) produces the Type 
pair in ([9]). 

Consider the Type I block AA„ with [A| = 1. If a matrix Fn is non- 
singular and F~*Fn is similar to Jn{^), then AA„ is *congruent to ±F„. 
Suppose A^ 7^ —1. Then 

_ . 1 - A^ _ . A(l - A^) _ . A- A _ ImA 
'~ 'l + A2 ~ 'a(1 + A2) ~ 'a + A ~ 

is real and 

^2 l + ic 
1 — tc 
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Consider the symmetric matrix 

Fn := A„(l + ic, i) = Pn {In + iJn{c)) 

Then 



(44) 



Fn 



{In - iJn{c)) ^ Pn 



and 



{In-iJn{c)) ^ PnPn{In + iJn{c)) 
= ((1 - ic) In - iJn{0)y^ ((1 + ic) In + iJn{0)) 

= A2 (/„ + ai J„(0) + a2 J„(0)2 + 03 J„(0)=^ + •••), 

in which oi = 2i {l + c^) ^ 7^ 0. Thus, F~*Fn is similar to Jn(A^), so AA„ is 
*congruent to =bA„(l + zc, i). The Cartesian decomposition of ibA„(l + ic, i) 
produces the first Type I pair in ([9]). 

Now suppose that = — 1 and consider the real matrix 



Gn '■- 







1 

-1 1 



-1 1 



m 



m 



if n = 2m, 



1 





1 



(45) 



m 



} 1 if n = 2m + 1 



m 



(the boxed unit is at the center). Since G~^Gn is similar to J„(— 1), AA„ 
is *congruent to ±G„. The Cartesian decomposition of ±G„ produces the 
second Type I pair in Q. 

The Type II block H2n(p) with \fi\ > 1 is *congruent to H2n{p'~^) be- 
cause their *cosquares are both similar to Jn(/i) ® Jn{p'~^)- Represent p.~^, 
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a point in the open unit disk with the origin is omitted, as 

1 l + iu 

1 — lU 

in which is in the open upper half plane with the point i omitted. In fact, 



2lm n + i i\fj,\ — 1 

z/ = ^ L —a + iby^i, b>0, a,beR. (46) 

l/i + 1[ 

We have 

(In + iJn{y)) {In " « JnM)~^ = ((1 + il^) In + Uni^)) ((1 " i^^) In " ? Jn(0))" 

= /i'l (/„ + ai J„(0) + 02 J„(0)2 + 03 J„(Of + • 
in which oi = 2i[l + u"^)'^ / 0. Thus, 

(/„ + iJnW)) {In - iJn{v))~^ 

is similar to Jn{p'~^)- Let be a nonsingular matrix such that 

S'"V„(/2"^)S' = (/„ + iJ„(i^)) {In - iJn{i^)T^ 
and compute the following *congruence of H2n{p'~^)- 



S{In-iJn{l^)) ■ 

S-* 



/„ 

Jn{fl'^) 



s(/„-iJ„H) 
s- 

In + iJn{l') 

In + iJn{v) 



The Cartesian decomposition of this matrix produces the Type II pair in 
([9]) with the parameters a and b defined in (j46p . □ 



For a given A with |A| = 1, the =b signs associated with the Type I canon- 
ical pairs in ([9]) can be determined using the algorithms in either Section 3 
or Section 4. 

For example, suppose n = 2 and X = i. The matrix 



Go 



1 
-1 1 
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defined in (|45p was analyzed in Example [H We found that iA2 is *congruent 
to — G2, so the *congruence canonical pair associated with iA2 is 





" 


■ 




" 


1 ■ 




( 





1 




1 





) 



(A2(0,1),A2(1,0)) 



As a second example, suppose n = 2 and |A| = 1, but / -1. Let 
A = a + z6 (a, 6 G M). The matrix 



F2 = A2(l + ib/a, i) = A2(A/a, i) 



A/a 
A/a i 



defined in (j44p . was analyzed in Example [9j We found that AA2 is *congru- 
ent to -F2 if a > 0, and to —F2 if a < 0. Thus, the *congruence canonical 
pair associated with AA2 is 



(A2(l,0), A2(6/a,l)) 
if Re A > 0, and is 

if Re A < 0. 





"01' 




" b/a' 




( 


1 




b/a 1 


) 





" 1 ■ 




" b/a' 


( 


1 




b/a 1 
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